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CLAIMS 94901 238 

A programmed computer' system for designing optimal 
oligonucleotide sequences for use with a gene sequence data 
source comprising : 

first input means for introducing user-selected gene 
sequence into the computer system ; 

memory means for storing user-selected ge ne sequence; 
means for accessing gene sequence data from -said gene 
sequence data source; 

means for performing exact and inexact match modelling 
between gene sequences; 

means for performing hybridization strength modelling on 
gene sequences; 

means for selecting either of said modelling means; and, 
means for presenting the results of said modelling to 

present candidate oligonucleotide sequences. — 

2~\ A programmed computer_ system ' i n accordance with claim 
1, wherein said means for performing exact and inexact match 
modelling utilizes said accessing means to introduce a user- 
selected set of gene sequence data and user-selected set of 
target gene sequence data from said gene sequence data source 
into the computer system and said memory means to store said 
gene sequence data and said target gene sequence data and 
wherein said meanis for performing exact and inexact match 
modelling includes: 

-means for determining a minimum sequence length; 
means for creating a look-up hash table and linked list 
in memory for each gene sequence in said gene sequence data 
and each of said garget gene sequences; 
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means for calculating the minimum length of any matching 
gene subsequence of said gene sequence data and said target 
gene sequence data; 

means for comparing each base pair character in each said 
target sequence stored in a hash table in memory to each base 
pair character of said gene" sequence stored in a hash table in 
memory; 

means for finding a matching seed by determining if the 
said comparison results in a matching gene subsequence of 
length equal to said calculated minimum length; 

means for comparing base pair characters behind and ahead 
of said seed to determine if there exists an extended match of 
a subsequence of base pair characters of length greater than 
the calculated minimum length, resulting in a current hit 
sequence ; 

means for calculating whether said current hit sequence 
is longer than said minimum sequence length, resulting in a 
current candidate oligonucleotide sequence; 

means for storing said current candidate oligonucleotide 

sequence ; and , 

wherein said presenting means provides said current 
candidate oligonucleotide sequence to the user. 

3,- A programmed computer system in accordance with claim 
2, wherein said computer system includes: 

means for calculating the melting temperature for each 
candidate oligonucleotide sequence; 

means for tracking the number and melting temperature of 
the matches for each candidate oligonucleotide sequence; 
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means for tracking the location of a set number of the 
best candidate oligonucleotide sequences preferably employing 
a priority queue by sorting said candidate oligonucleotide 
sequence in reverse order and sorting said oligonucleotides by 
hybridization strength; and, 

wherein said presenting means is operative to present 
said additional results to the user; and, 

wherein said presenting means operably provides said 
melting temperature to the user. 

4. A programmed computer system in accordance with 
claims 2 or 3, wherein said first input means in operative to 
introduce a user-selected maximum number of mismatches and a 
user-selected minimum candidate oligonucleotide sequence 
length into the computer system, and wherein said means for 
calculating the minimum length of any matching gene 
subsequence of said gene sequence data and said target gene 
sequence data comprises the steps of: 

means for subtracting said maximum number of mismatches 
from said minimum candidate oligonucleotide sequence length to 

give a first result; 

means for dividing said first result by said maximum 
number of mismatches plus one to give a second result; 

means for incrementing said second result by one if the 
remainder is not equal to zero to give a third result; and, 

means for truncating said third result to an integer, and 
preferably said means for calculating the hairpin 
characteristics of said candidate oligonucleotide sequence 
comprises the steps of: 
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calculating a complenientary sequence to the candidate 
oligonucleotide sequence by reversing the base pair order of 
the candidate oligonucleotide and substituting complementary 
base pairs; 

comparing each character of said original candidate 
oligonucleotide sequence and said complementary sequence; 

finding the longest match between said original candidate 
oligonucleotide sequence and said complementary sequence; and, 

saving the match with the longest hairpin distance if any 
two matches have the same length; 

means for storing hairpin characteristics; and, 

wherein said presenting means provides said hairpin 
characteristics to the user. 

5. A programmed computer system in accordance with any 
preceding claim, wherein said means for performing 
hybridization strength modelling utilizes said first input 
means to introduce a user -se 1 ect ed screening threshold into 
the computer system and said accessing means to introduce a 
user-selected set of gene sequence data from said gene 
sequence data source into the computer system, and said memory 
means to store said gene sequence data, said target gene 
sequence data and said screening threshold and wherein said 
means for performing hybridization strength modelling 
comprises : 

means for preprocessing said target gene sequence data 
and said gene sequence data by selecting only those sequences 

without introns; 

means for forming a preparation file of gene sequence 
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fragments by cutting said target gene sequences into fixed 
length target gone subsequences and sorting said subsequences 

in lexicographical order; 

means for merge sorting said gene sequences; 

means for forming multiple lists of screens by forming 
lists of subsequences of the preparation file of length equal 
to said screening threshold; 

means for indexing, sorting and storing said screens in 

said memory means; 

means for sequentially comparing said preparation file 
gene sequences with each of said screens to design candidate 
oligonucleotide sequencies; 

means for calculating the hybridization strengths between 
a gene sequence and all candidate oligonucleotide sequences 
containing that gene sequence by accounting for Guanine- 
Cytosine (GC) and Adenine-Thymine (AT) base pair content of 
the gene sequence and the number of mismatches between said 
preparation file sequences and a said screen when said 
comparison results in a match; 

means for preparing the candidate oligonucleotide 
sequence and hybridization strength for presentation to the 
user ; and 

wherein said presenting means provides the candidate 
oligonucleotide sequence and hybridization strength to the 
user- 

6, A programmed computer system in accordance with claim 
5, wherein said computer system includes: 

means for assigning a numerical score to each said gene 
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sequence; and 

means for sorting said gene sequences in accordance with 
said numerical score, and/or preferably wherein means for 
assigning a numerical score to each said gene sequence 
operates by tallying the quantity "exp" where "exp" « > e"*^ 

and wherein Tm is the melting temperature for the said gene 
sequence; and 

means for sorting said gene sequences in accordance with 
said numerical score. 

7. A programmed computer system in accordance with claim 
5, wherein said means for calculating the hybridization 
strengths between a gene sequence and all candidate 
oligonucleotide sequences containing that gene sequence 
comprises the steps of: 

accessing gene sequence data from said gene sequence data 

source; 

comparing base pairs of a first gene sequence and a 
second gene sequence to determine if <^match existsT) 

incrementing sai d first gene se gu^ce^s^bo^^^^ strength by^ 
some first numberif a base pai r _cha£act er in said ^ 
sequence and said second gene sequence matc h and the matched 
t^^e^'parr is equal to a combination of the bases Guanine [G) 
and Cytosine (C); 

incrementing said first gene sequence's bound strength by 
some second number if a base pair character in said first gene 
sequence and said second gene sequence match and the matched 
base pair is equal to a combination of the bases Adenine (A) 
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and Thymine (T) ; 

decrementing said first gene sequence's bound strength by 
a third number if there is no match in base pairs between said 
first gene sequence and said second gene sequence; 

comparing said first gene sequence's bound strength to 
said first gene sequence's unbound strength; 

setting said first gene sequence's unbound strength equal 
to its bound strength if said first gene sequence's bound 
strength is greater than said first gene sequence's unbound 
strength; and 

resetting said first gene sequence's bound strength to 
zero if said first gene sequence's unbound strength is less 
than zero- 

8. A programmed computer system in accordance with claim 
5, wherein said computer system includes a means for 
calculating the hairpin characteristics of said candidate 
ol igonucleot ide sequence; 

means for preparing the hairpin characteristics for 
presentation to the user; and, 

wherein said presenting means provides the hairpin 
characteristics to the user, wherein preferably means for 
calculating the hair pin characteristics of said candidate 
oligonucleotide sequence comprises the steps of: 

calculating a complementary sequence to the candidate 
oligonucleotide sequence by reversing the base pair order of 
the candidate oligonucleotide sequence and substituting 
complementary base pairs; 

comparing each character of said original candidate 
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oligonucleotide sequence and said complementary sequencer- 
finding the longest match between said original candidate 
oligonucleotide sequence and said complementary sequence; and 

saving the match with the longest hairpin distance if any 
two matches have the same length; 

means for preparing the hairpin characteristics for 
presentation to the user; and 

wherein said presenting means provides the hairpin 
characteristics to the user. 

9, A programmed computer system in accordance with claim 
1, wherein said means for performing exact and inexact match 
modelling utilizes said accessing means to introduce a user- 
selected set of gene sequence data and user -se 1 ect ed set of 
target gene sequence data from said gene sequence data source 
into the computer system and said memory means to store said 
gene sequence data and said target gene, sequence data and 
wherein said means for performing exact and inexact match 
modelling includes: 

means for determining a minimum sequence length; 
means for creating a look-up hash table and linked list 
in memory for each gene sequence in said gene sequence data 
and each of said target gene sequences; 

means for calculating the minimum length of any matching 
gene subsequence of said gene sequence data and said target 
gene sequence data; 

means for transforming base characters in each said 
target sequence and in each said gene sequence into numeric 
digits; 
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means for comparing each base pair digit in each said 
target sequence stored in a hash table in memory to each base 
pair digit of said gene sequence stored in a hash table in 
memory; 

means for finding a matching seed by determining if the 
said comparison results in a matching gene subsequence of - 
length equal to said calculated minimum length; 

means for comparing base pair digits behind and ahead of 
said seed to determine if there exists an extended match of a 
subsequence of base pair digits of length greater than the 
calculated minimum length, resulting in a current hit 
sequence; 

means for calculat ing whether ^said current hit sequence 
is longer than said minimum sequence length, resulting in a 
current candidate oligonucleotide sequence; 

means for storing said current candidate oligonucleotide 

sequence; and 

wherein said presenting means provides said current 
candidate oligonucleotide sequence to the user. 

A programmed computer system for designing candidate 
olig^rae^ieotide sequences for use with a gene sequence data 
source including: 

first input means for introducing user-selected gene 
sequence, design, model and presentation criteria and user- 
specified sequence length into the computer system; 

memory means for storing said gene sequence, design, 
model and presentation criteria and said sequence length; 

means for accessing gene sequence data from said gene 
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sequence data source; 

wherein said accessing means is operative to introduce a 
user-selected set of gene sequence data and a user-selected 
set of target gene sequence data from said gene sequence data 
source into the computer system; 

wherein said criteria are used for comparison of gene 
sequence data and target gene sequence data; 

means for comparing said gene sequences against said 
target gene sequences employing said criteria; 

means for calculating candidate oligonucleotide sequences 
of said sequence length that are either common to a pool of 
user-specified gene sequences or specific to a particular 
user-specified gene sequence; 

means for calculating the homology between the candidate 
oligonucleotide sequences and said gene sequence data; 

means for calculating a candidate oligonucleotide 
sequence's hairpin characteristics; 

means for displaying in multiple dimensions the gene 
sequences which result from the comparisons and calculations 
characterized in that said display format exhibits: 

the starting position of each candidate oligonucleotide 
sequence in one dimension; 

a candidate oligonucleotide sequence's specificity to the 
target gene sequence in a second dimension; and 

superimposed melting temperatures of gene sequences in 
contrasting presentations in at least an apparent third 
dimension; 

wherein said display further includes a cursor movable 



along one dimension of said display that selects a position 
for an expansion of data representing the homology between the 
candidate oligonucleotide sequences and said gene sequence 
data; and 

wherein said display is operative to display in 
alphanumeric form the homology between the candidate 
oligonucleotide sequences and said gene sequence data; and 

wherein said display is operative to provide an expansion 
of data including presenting 

false hybridizations at various melting temperatures for 
all candidate oligonucleotide sequences; 

the location of each false hybridization; 

a candidate oligonucleotide sequence's starting position; 

and, 

hairpin characteristics of each candidate oligonucleotide 
siequence. 

11. A method for designing candidate oligonucleotide 
sequences by performing exact and inexact match modelling for. 
use with a gene sequence data source comprising the steps of: 

introducing user-selected gene sequence into a computer 
system; 

accessing gene sequence data from said gene sequence data 
source ; 

storing user-selected gene sequence in the memory of the 
computer system; 

accessing the gene sequence source to introduce the user- 

r " 

selected set of gene sequence data and user-selected set of 
target gene sequence data from said gene sequence data source 
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into the computer system; 

storing said gene sequence data and said target gene 
sequence data in the memory of the computer system; 

determining a minimum sequence length; 

creating a look-up hash table and linked list in memory 
for each gene sequence in said gene sequence data and each of 
said target gene sequences; 

calculating the minimum length of any matching gene 
subsequence of said gene sequence data and said target gene 
sequence data; 

comparing each base pair character in each said target 
sequence stored in a hash table in memory to each base pair 
character of said gene sequence stored in a hash table in 
memory; . 

determining a matching seed by determining if the said 
comparison results in a matching gene subsequence of length 
equal to said calculated minimum length; 

comparing base pair characters behind and ahead of said 
seed to determine if there exists an extended match of a 
subsequence of base pair characters of length greater than the 
calculated minimum length, resulting in a current hit 
sequence ; 

calculating whether said current hit sequence is longer 
than said minimum sequence length, resulting in a current 
candidate ol igonucleot ide sequence; 

storing said current candidate oligonucleotide sequence 
in the memory of the computer system; and, 

presenting a representation of said current candidate 
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oligonucleotide sequence to the user. 

12. A method in accordance with claim 11, wherein said 
method includes the steps for performing additional 
calculations for each candidate oligonucleotide probe, said 
additional calculations comprising: 

calculating -the melting temperature for each candidate 
oligonucleotide sequence; 

tracking the number and melting temperature of the 
matches for each candidate oligonucleotide sequence; 

tracking the location of a set number of the best 
candidate oligonucleotide sequences preferably employing a 
priority queue by sorting said candidate oligonucleotide 
sequences in reverse order and sorting said candidate 
oligonucleotide sequences by hybridization strength; and 
presenting said additional results to the user. 
13. A method in accordance with claim 12, wherein said 
step for calculating the minimum length of any matching gene 
subsequence comprises: 

introducing a user-selected maximum number of mismatches 
and a user-selected minimum candidate oligonucleotide sequence 
length into the computer system; 

subtracting said maximum number of mismatches from said 
minimum candidate oligonucleotide sequence length to give a 
first result; 

dividing said first result by said maximum number of 
mismatches plus one to give a second result; 

incrementing said second result by one if the remainder 
is not equal to zero to give a third result; and. 
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truncating said third result to an integer. 
14. A method in accordance with any one of claims 11 to 
13, wherein said method includes the step of calculating the 
hairpin characteristics of said candidate oligonucleotide 
sequence, preferably said method includes the step' of 
calculating the hairpin characteristics of said candidate 
oligonucleotide sequence comprising: 

calculating a complementary sequence to the candidate 
oligonucleotide sequence by reversing the base pair order of 
the candidate oligonucleotide sequence and substituting 
complementary base pairs; 

comparing each character of said original candidate 
oligonucleotide sequence and said complementary sequence; 

finding the longest match between said original candidate 
oligonucleotide sequence and said complementary sequence; and, 
saving the match with the' longest hairpin distance if any 
two matches have the same length, , 

/ui5,j))a method for designing candidate oligonucleotide 
sequences by performing hybridization strength modelling for 
use with a gene sequence data source comprising the steps of: 

introducing user-selected gene sequence and a user- 
selected screening threshold into a computer system; 

storing user-selected gene sequence and said screening 
threshold in the memory of the computer system; 

accessing the gene sequence source to introduce the user- 
selected set of gene sequence data and a user-selected set of 
target gene sequence data from said gene sequence data source 
into the computer system; . 
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storing said gene sequence data and said target gene 
sequence data in the memory of the computer system; 

preprocessing said target gene sequence data and said 
gene sequence data by selecting only those sequences without 
introns ; 

forming a preparation file of gene sequence fragments by 
cutting said target gene sequences into fixed length target 
gene subsequences and sorting said subsequences in 
lexicographical order; 

merge sorting said gene sequences; 

forming multiple lists of screens by forming lists of 
subsequences of the preparation fiie of length equal to said 
screening threshold; 

indexing and sorting said screens in memory; 

storing said screens in the memory of the computer 
system ; 

sequentially comparing said preparation file gene 
sequences with each of said screens to design candidate 
oligonucleotide sequences; 

calculating the hybridization strengths between a gene 
sequence and all candidate oligonucleotide sequences 
containing that gene sequence by accounting for Guanine- 
Cytosine (GC) and Adenine -Thymine (AT) base pair content of 
the gene sequence and the number of mismatches between said 
preparation file sequences and a said screen when said 
comparison results in a match; 

preparing the candidate oligonucleotide sequence and 
hybridization strength for presentation to the user; and. 
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presenting the candidate oligonucleotide sequence and 
hybridization strength to the user, and preferably wherein 
said method includes the steps for performing additional 
calculations for each candidate oligonucleotide probe, said 
additional calculations comprising: 

calculating the melting temperature for each candidate 
oligonucleotide sequence; 

tracking the number and melting temperature for each 
candidate oligonucleotide sequence; 

tracking the location of a set number of the best 
candidate oligonucleotide sequences preferably employing a 
priority queue by sorting said candidate oligonucleotide 
sequences in reverse order and sorting said candidate 
oligonucleotide sequences by hybridization strength; and 

presenting said additional results to the user. 

16. A method in accordance with claim 15, to use with a 
gene sequence data source, programmed to determine 
hybridization strength comprising the steps of: 

comparing base pairs of a first gene sequence and a 
second gene sequence to determine if a match exists; 

incrementing said first gene sequence's bound strength by 
some first number if a base pair character in said first gene 
sequence and said second gene sequence match and the matched 
base pair is equal to a combination of the bases Guanine (G) 

and Cytosine (C); 

incrementing said first gene sequence's bound strength by 
some second number if a base pair character in said first gene 
sequence and said second gene sequence match and the matched 
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base pair is equal to a combination of the ba.es Adenine (A) 

and Thymine (T) ; 

decrementing said first gene sequence's bound strength by 
a third number if there is no match in base pairs between said 
first gene sequence and said second . gene ..^eqiience; 

comparing said first gene sequence's bound strength to 
said first gene sequence's unbound strength; 

setting said first gene sequence's unbound strength equal 
to its bound strength if said first gene sequence's bound 
strength is greater than said first gene sequence's unbound 
strength; and 

resetting said first gene sequence's bound strength to 
^ero if said first gene sequence's unbound strength is less 

than zero; and 

preferably wherein said first and second numbers are 

greater than zero; and/or 

preferably wherein said second number is in the order of 

42* of said first number; and/or 

wherein said third number is in the order of 5% larger 
than said first number; and/or 

preferably wherein said method includes the step of 
calculating the hairpin characteristics of said candidate 
oligonucleotide sequence, wherein preferably including the 
steps of: 

calculating a complementary sequence to the candidate 
oligonucleotide sequence by reversing the base pair order of 
the candidate oligonucleotide sequence and substituting 
complementary base pairs; 
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comparing each character of said original candidate 

oligonucleotide sequence and said complementary sequencer- 
finding the longest match between said original candidate 

oligonucleotide sequence and said complementary sequence; and, 
saving the match with the longest hairpin distance if any 

two matches have the same length. 

(vlP^ A method for designing candidate oligonucleotide 

sequences for use with a gene sequence data source comprising 

the steps of: 

introducing user-selected gene sequence and user- 
specified sequence length into a computer system; 

storing said gene sequence and said sequence length in 
the memory of the computer system; 

accessing gene sequence data from said gene sequence data 
source; 

accessing the gene sequence source to introduce the user- 
selected set of gene sequence data and a user-selected set of 
target gene sequence data from said gene sequence data source 
into the computer system; 

comparing said gene sequences against said target gene 
sequences employing said cr i t er i a ; 

calculating candidate oligonucleotide sequences of said 
sequence length that are either common to a pool of user- 
specified gene sequences or specific to a particular user- 
specified gene sequence; 

calculating the homology between the candidate 
oligonucleotide sequences and said gene sequence, data; 

displaying in multiple dimensions the gene sequences 
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which result from the comparisons and calculations 
characterized in that said display format exhibits; 

the starting position of each candidate oligonucleotide 

sequence in one dimension; 

a candidate oligonucleotide sequence's specificity to the 
target gene sequence in a second dimension; and, 

superimposed melting temperatures of gene sequences in 
contrasting presentations in at least an apparent third 
dimension, and wherein preferably said method includes the 
step of calculating a candidate oligonucleotide sequence's 
hairpin characteristics, preferably comprising the steps of: 

calculating a complementary sequence to the said gene 
: sequence by reversing the base pair order of the gene sequence 
and substituting complementary base pairs; 

comparing each character of said original gene sequence 
and said complementary sequence; 

finding the longest match between said original gene 
sequence and said complementary sequence; and, 

saving the match with the longest hairpin distance if any 
two matches have the same length - 

18. A method to determine hybridization strength between 
two or more gene sequences for use with a gene sequence data 
source, comprising the steps of: 

accessing gene sequence data from said gene sequence data 

source; . 

comparing base pairs of a first gene sequence and a 
second gene sequence to determine if a match exists; 

incrementing said first gene sequence ' s bound strength by 
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some first number if a base pair character in said first gene 
sequence and said second gene sequence match and the matched 
base pair is equal to a combination of the bases Guanine (G) 
and Cytosine (C) ; 

incrementing said first gene sequence's bound strength by 
some second number if a base pair character in said first gene 
sequence and said second gene sequence match and the matched 
base pair is equal to a combination of the bases Adenine (A) 

and Thymine (T) ; 

decrementing said first gene sequence's bound strength by 
a third number if there is no match in base pairs between said 
first gene sequence and second gene sequence; 

comparing said, first gene sequence's bound strength to 
said first gene sequence's unbound strength; 

setting said first gene sequence's unbound strength equal 
to its bound strength if said first gene sequence's bound 
strength is greater than said first gene sequence's unbound 
strength; and 

resetting said first gene sequence's bound strength to 
zero if said first gene sequence's unbound strength is less 
than zero, and wherein preferably said first and second 
numbers are greater than zero, and/or preferably wherein said 
second number is in the order of 42* of said first number, 
and/or wherein said third number is in the order of 5% larger 
than said first number. 

19. A method of calculating the minimum length of any 
matching gene subsequence comprising: 

introducing a user-selected maximum number of mismatches 
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and a user-selected minimum candidate, ol igonucleot ide sequence 
length; 

subtracting said maximum number of mismatches from said 
minimum candidate oligonucleotide sequence length to give a 
first result ; 

dividing said first result by said maximum number of 
mismatches plus one to give a second result; 

incrementing said second result by one if the remainder 
is not equal to zero to give a third result; and, 

truncating said third result to an integer. 

20. A method of calculating ^^airpin c:haracterj_s^^^ 
a gene sequence comprising: 

calculating a complementary sequence to the said gene 
sequence by reversing the base pair order of the gene sequence 
and substituting complementary base pairs; 

comparing each character of said original gene sequence 
and said complementary sequence; 

finding the longest match between said original gene 
sequence and said complementary sequence; and, 

saving the match with the longest hairpin distance if any 
two matches have the same length. 



